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ABSTRACT 

Previous calculations of the pregalactic chemistry have found that a small amount 
of H2, x[H2] = n[H2]/n[H] « 2.6 x 10~^, is produced catalytically through the H~, 
H^, and HeH"*" mechanisms. We revisit this standard calculation taking into account 
the effects of the nonthermal radiation background produced by cosmic hydrogen 
recombination, which is particularly effective at destroying H~ via photodetachment. 
We also take into consideration the non-equilibrium level populations of H^, which 
occur since transitions among the rotational-vibrational levels are slow compared to 
photodissociation. The new calculation predicts a final H2 abundance of x[H2] ~ 
6 X 10^^ for the standard cosmology. This production is due almost entirely to the 
H~ mechanism, with ~ 1 per cent coming from HeH+ and ~ 0.004 per cent from H^. 
We evaluate the heating of the diffuse pregalactic gas from the chemical reactions that 
produce H2 and from rotational transitions in H2, and find them to be negligible. 

Key words: cosmology: theory - intergalactic medium - molecular processes. 



1 INTRODUCTION 

, One of the key problems in cosmology is to understand 
■ the physical and chemical state of the baryonic matter in 
' the Universe. At high redshift, the baryonic matter was 
, fully ionized and co-existed with a thermalized radiation 
' field (the cosmic microwave background, or CMB). By red- 
shift z ~ 10^, the Universe had expanded and cooled to 
~ 3000 K, at which point the ionized nuclei and free elec- 
trons of the primordial plasma combined to form neutral 
atoms. T his cosm i c reco mbina tion was first stu d ied th eoret- 
ically by iPeebled lll968D and IZel'dovich eTsIl Jl968D . The 
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combination happened, and that it occurred over a narrow 
range in redshift, in accordance with predictions. 

As the Universe continued to expand and cool, the for- 
mation of molecules became thermodynamically favourable. 
Since hydrogen is most abundant, one would expect the 
most abundant molecule to be H2. However, unlike atomic 
recombination, which occurs shortly after it becomes ther- 
modynamica lly favourable and proceeds nearly to com- 
pletion (e.g. ISeaeer et al.ll2000l) . cosmological formation of 
molecules is slow and freezes out with a final abundance, 
a;[H2] = n[H2]/n[H] < 1. 



Despite their small abundance, molecules in the early 
universe have been investigated for several reasons. The first 
is that the primordial gas, mainly of hydrogen and helium 
atoms, lacks the low-lying excitations necessary for cool- 
ing and therefore star formation at low temperatures. On 
the other hand, molecules (which possess low-lying rota- 
tional excitations) could provide the coo ling necessary to 
form the first stars dSaslaw fc Zipovlll96i1) . However, recent 
calculations indicate that the primordial H2 abundance is 
far too small for this, and that the only H2 important for 
coohng of early hal oes is formed in collapsed haloes (e.g. 
iTegmark et al.lll997l) . A second reason for studying H2 pro- 
duction is that the heating of the gas, either via rotational 
transitions induced by the CMB or the chemical energy 
released by formation of the molecules, could affect the 
temperature of the pregalactic gas. Here even a small ef- 
fect could be important for proposals to study the absorp- 
tion of the CMB by prega lactic gas in the H I 21-cm line 
jLoeb fc Zaldarriagall20M) . Finally, there is the (perhaps 
academic) motivation to understand the composition of the 
primordial gas as part of elucidating the standard cosmolog- 
ical model. 

Th e first calculatio n of th e primordial H2 abundance 
was by ISaslaw fc Zipovl I^G^)- Noting that the direct ra- 
diative association of two H atoms is forbidden, they pro- 
posed that H2 molecules could be built up using HJ as 
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an intermedia te state (the specific re actio ns will be Riven 
in Sec tion I^J . iPeebles fc Dickj il968l) and iHirasawa et alJ 
il969l) suggested that the H~ mechanism dominated the 
production of molecules in primordial gas clouds. A num- 
ber of subsequent studies considered in increasing detail the 
H2 abundance a,nd cooling in primordial clouds ifH irasawal 
ll969l:lYonevamdll972l:lHutchinsll97^ . lLepp fc ShulJ lll984l) 
performed a calculation of the abundances of H2 as well as 
HD and LiH, calculating a final abundance a;[H2] ~ 10~^ in 
the intergalactic gas. This is essentially today's "standard" 
calculation of the primordial H2 abundance, although some 
of the reaction rates and cosmological parameters have been 
updated. These updated analyses, which include substantial 
revision s to the deuterium and lithium chemistry, can be 
found i nlPuv et alJ ^1993^■lPalla et alJ lll995H . lGalli fc Fallal 
il998h . l^£mjire^]T ' i^998f) . and references therein. 



The effect of H2 on h eating of t h e gas before collapse 
has also been considered. IPuv et al.l lll993i) presented the 
first analysis of the thermal effect of molecules; they found 
a moderate effect (~ 10 per cent) on the gas tempera- 
ture mainly due to rotational lines in H2, and a smaller 
effect due to chemical reactions. The heating i n rota- 
tion al lines was al s o con sidered by iHaiman et alJ (1199681) 
and iGalli fc Pallal (|l99^, b ut was revised downward by 
iFlower &: Pineau des Foretj i200Cf) . who concluded the ef- 
fect was insignificant. We revisit the chemical heating here, 
and conclude that it dominates over H2 rotational lines. 
However, even this effect is probably too small to be de- 
tected by 21-cm experiments in the forseeable future as it 
changes the gas temperature at the ~ 10"* level. 

All of these analyses, however, have been based 
on several common assumptions. One is the assumption 
of a purely thermal radiation field, which is not com- 
pletely correct because of the line and continuum radiation 
emitted during hydrogen and helium recombina tion (e.g. 
iRvbicki & Dell'Antonio 'l99j; IWong et all l2006f) . Indeed, 
[Switzor fc _Hira ta (2005) found that this spectral distortion 
suppresses the lithium abundance by several o rders of mag- 
nitude as compared with previous calculation s iStancil et alJ 
ll996l:lGalh fc Pallalll993lStancil et al.ll2002h . The other as- 
sumption is the use of photodissociation rates based ei- 
ther on local thermodynamic equilibrium (LTE) populations 
of the rotational-vibrational levels, or with all ions in 
the ground state. It is however known that forms pref- 
erentially in excited states and since radiative transitions 
between levels are slow, LTE m ay not apply. The im por- 
tance of th i s was recognized by iGalli fc Pallal (|^9^ and 
iLepp et alJ i2002l) . but a full non-LTE analysis of level 
populations has not been done. 

Our purpose here is to revisit the calculation of H2 
abundance, including spectral distortions to the CMB and 
with a level-resolved treatment of HJ. In Section |51 we in- 
troduce the chemical reactions important for H2 production. 
The H~ mechanism, including the effect of the spectral dis- 
tortion, is discussed in Section |21 The H J mechanism and 
the level-resolved treatment is in Section |3] and HeH"*" is 
discussed in Section]^ The abundances of H2, H~, Hj, and 
HeH"'' are calculated for our presently favoured cosmology 
in Section 1^ The heating of the pregalactic gas by H2 and 
the chemical reactions leading to its formation is considered 
in Section |7| We conclude in Section |H| The theory of the 



energy levels and transitions of the ion is recapitulated 
in Appendix IXI 

In this paper, we have assumed a primordial helium 
abundance of Yp = 0.24, and a fl at AC DM cosmology 
with parameters from ISeliak et alJ (l2005l) : Q.b = 0.0462, 

= 0.281, and Ho = 71.0kms"^ Mpc"^ The number 
density n will refer to the total proper density of hydro- 
gen nuclei in all forms (ionized, atomic, and molecular), al- 
though in the regime of interest here it is mostly atomic. 
The notation xt (or x[i] for H~, Hj, and HeH^) will denote 
the number density of species i relative to the total num- 
ber density of hydrogen nuclei in all chemical forms (e.g. 
k[H2] = n[H2]/n = 1/2 if all hydrogen is molecular). 



2 THE REACTIONS 

Due to the lack of a dipole moment, it is forbidden for two 
H atoms to combine radiatively to form H2. Therefore cos- 
mological H2 production proceeds thr ough two main mecha- 
nisms catalyzed by cha,rKed particles (ISaslaw fc Zipovlll967l : 
iPeebles fc Dickelll968t lHirasawalll969r) . The H~ mechanism 
begins with radiative attachment to form H~, 

H + e-^H-+7, (1) 

and is completed when the associative detachment reaction 

H-+H^H2+e-. (2) 

A minor reaction that can use up H~ ions is mutual neu- 
tralization. 



2H. 



(3) 



An alternative mechanism is via Hj, in which the cat- 
alyst is a proton rather than an electron. 



H + H^ 



H++7. 



(4) 



The ion can then be converted to H2 via the reaction 



H ^ H2 



(5) 



It is also possible for HJT to be destroyed by dissociative 
recombination, 



2H. 



(6) 



A third route is via the HeH mechanism. This is urnm- 
portant in the "standard" calculation, but given that we are 
revising the H~ and Hji" rates downward it is only prudent 
to include it. It begins with the production of HeH^ by ra- 
diative association, 



He H+ ^ HeH+ + 7. 



(7) 



This ion could either be photodissociated (reverse of Eq.|7|l, 
but it could also form Hj by the reaction 

HeH+ -hH^ He-hHj. (8) 

The HJ ion then participates in the usual sequence of reac- 
tions, Eqs. (gHHJ. 

All three of these mechanisms are suppressed at high 
redshift due to the intense CMB radiation, which drives 
Eqs. (0, and ^ strongly to the left. They are also 
inefficient at low redshift because the collisions required for 
them to proceed become rare as the universe expands and 
the gas density drops. Past calculations have H2 production 
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peaking at z ~ 260 for the Hj" mechanis m and 2; ~ 90 for 
the H~ mechanism (e.g. iLepp et al.l|2002h . 

Once formed, the H2 molecule can be destroyed by UV 
phot odissociation ( Stcchcr & Williams' 'l967l: lAberall et alJ 
ll992 : .Haiman et al. 1996b). Photooxcitation of H2 from the 
ground X^'E'^ electronic state to either B^Sj (Lyman band 
absorption) or C^llu (Werner band absorption), is usually 
followed by radiative decay back to the X^T,'^ electronic 
state. However, it is possible that this process leaves the 
H2 molecule in an unbound vibrational state, resulting in 
its dissociation into two H atoms (a table of prob abili- 
ties can be found in e.g. iDakarno fc Stephens! [l97cl) . Of 
course, this relies on the existence of UV radiation in the 
Lyman and Werner bands, which have minimum energies of 
11.2 and 12.3 eV respectively. During the postrecombination 
era, these energies are 2-3 orders of magnitude larger than 
kTcMB ~ 0.2 eV, and so the photodissociation of H2 plays a 
negligible role in determining the pregalactic H2 abundance. 
The presence of the spectral distortion does not change this 
situation since the distortion extends only up to 10.2 eV 
(the H I Lya energy). It is only after the first astrophysical 
sources turn on that intergalactic II2 can be destroyed. 

The cosmological production of H2 can be followed by 
keeping track of the abundances of the relevant species: H, 
He, H+, e", H", H+, HeH+, and Ha. The evolution of H, 
He, H^, and e~ has been investigated in the context of the 
cosmic recombination, and is essentially unaffected by the 
rates of the catalytic reactions Eqs. due to the small 

amount of H2 produ ced. We therefore use th e recombina- 
tion code Recfast ijSeager et alJ I1999I . l200(t to compute 
the abundances of H, H"*", and e~. The rate equation for H2 
is 



d3:[H2 
dt 



S-+S+, 



(9) 



where S- and S+ are the formation rates via the H~ and 
combined H^/HeH"'' mechanisms, respectively (note that 
the latter two cannot truly be separated since Eq. |H| cou- 
ples them). In the next several sections, we compute the H2 
production and photodissociation rates. Section|21treats the 
computation of S- (the H~) mechanims. Section 2] treats 
the computation of S+ including only the H J mechanism for 
simplicity. In Section|S]we introduce HeH"*" into the compu- 
tation of S+. 



3 THE H" MECHANISM 

The rate for the H~ mechanism depends on the population 
of H~ as a function of time. The relevant rate equation is 



da;[H- 
dt 



— kiXeXmn — fc_ia;[H ] — fc2a;[H ]a;Hifi 
-fc3a;[H"]a;[H+]n, 



(10) 

where fei is the rate for Eq. is the rate for the 

reverse reaction, k2 is the rate for Eq. Q, and fcs is the 
rate for Eq. l El. Th e forward rates are given by the fits in 
IStancil et aP lll998D : 

fci(r™) = 3xl0-^«(|^)"'%-^'"/«^^%m^-\ 

. rp . -0.1 

kiiTm) = l-SxlO'^f^j cm^s-\ and 



k3{T„ 



= 4 X 10" 



Vsooy 



(11) 



where the matter temperature, Tm, is in Kelvin. (The stim- 
ulated radiative attachment can be neglected in compari- 
son with the spontaneous rate because the energy of the 
emitted photon is always at least the H~ binding energy, 
or > 0.754 eV; this is much greater than fcsTcA/s in the 
redshift range of interest.) There is some uncertainty in the 
associative detachment rate k2 and the mutual neutraliza- 
tion rate fca, which we will discuss in Section]^ The ratio 
of H~ destruction rate to the Hubble rate is always at least 
fc2n/7i" « 3500[(l-f 2;)/10]^'^ > 1, so we may treat a; [H"] by 
the steady-state approximation, 



kiXeXmn 



k2xmn + fc-i + fcsxfH^ 



The production rate for H2 via H is then, 
7 7 22 

KiK2XeX^in 



S. 



k2Xmn + k-i + k3x[H+]n 



(12) 



(13) 



The photodetachment rate fc_i depends on the details 
of the radiation field and can be broken into thermal (black- 
body CMB) and nonthermal (spectral distortion) parts: 



fcL'?'(TcAis) 



(14) 



The rate from the thermal photons can be computed via the 
principle of detailed balance. 



k^'^\TcMB) 



rriekBTcMB 
y<ki{TcMB), 



3/2 



-B(H-)/kBTcMI 



(15) 



where the factor of 4 comes from the spin degeneracy (1 for 
H", 2 each for H and e"), and B(H") = 0.754 eV is the 
photodetachment threshold energy. Note that Tcmb is used 
here instead of Tm since photodetachment depends only on 
the properties of the radiation field. 

In the standard calculation the thermal rates are used, 
i.e. k-i = k^2i ^ ■ The nonthermal contribution to the pho- 
todetachment rate is 



B(U-)/h 



V 



(16) 



where r{v) is the number of distortion photons per H atom 
per logarithmic range in frequency, and (j-\{v) is the pho- 
todetachment c ross section, for which we use the fit by 
iTegmark etHI Jl99'i) . 



a = 3.486 X 10~ 



; {x - 



N3/2 



(17) 



where x = hv/B{¥{ ). The variable r{u) is related to the 
phase space density fiy) of photons by 



/(^) 



nc , , 



(18) 



where the first term denotes the thermal CMB contribution, 
and the latter term is the spectral distortion. The spectral 
distortion dominates at high frequencies, hv /kBTcMB > 30. 

The spectra l dist ortion r{i>) is calculated as in 
ISwitzer fc Hiratal lj200!t : the H I 2s ^ Is two-photon de- 
cay and Lyman-a resonance escape rates were obtained 
from Recfast and integrated as described in Section II 
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of ISwitz er fc Hirata| l|20 05l). We note that a recent compu- 
tation by IWong et alJ i20Ml find an additional distortion 
due to He l 2^^5*0 l^So two-photon decays and escape 
from the 2^P° — » 1 ^5'o resonance. As can be seen from 
Fig. 3 of IWong et all J^JO&'l . the He I distortion contributes 
significantly to the photon spectrum at short wavelengths 
A < 140(1 + z)~^ fim, and could be a significant contribu- 
tion to H~ or HJ destruction rates. However, these photons 
are in the H I Lyman continuum (A < 912 A) at z ~ 1500, 
when the universe is already optically thick (r ~ 10^) due to 
H I photoionization. Thus we do not expect the high-energy 
spectral distortion from He I recombination to survive and 
have not included it in our analysis. 



4 THE H+ MECHANISM 

The HJ mechanism is more complicated to analyze than 
the H~ mechanism because unlike H~, the H^ ion has many 
bound states. One must therefore determine the populations 
of each Hj level, taking into account the radiative associa- 
tion rates to each level, the radiative and coUisional rates for 
changing the rotational and vibrational quantum numbers, 
and the destruction rates by photodissociation and charge 
transfer. We start with an overview of the physics of the HJ 
ion, and describe the full reaction network that describes 
the ion. Finally, we show how the relevant physics can be 
captured by a "two-level" approximation; this provides a 
computationally simpler approach to the Hj mechanism. 

The HJ ion is diatomic, and can be described by spec- 
ifying the electronic state, a rotational quantum number 
describing the total orbital angular momentum, and a vi- 
brational quantum number v equal to the number of radial 
nodes. The electronic states of interest to us are the ground 
state IsCTg (X^S^) and the first excited state 2pau {A^T,i); 
higher states are not accessible at the temperatures under 
consideration. We ignore the spin-orbit coupling and hyper- 
fine structure since their energy splittings are small com- 
pared to fesTm or kBTcMB- Thus the degeneracy of a given 
rotational-vibrational level is g^N = 2{2N + l)(/nuc) where 
the 2 comes from the electron spin. The nuclear degeneracy 
is forced by proton wave function antisymmetry considera- 
tions to be = 1/4 for spatially symmetric states and 
ffnuc = 3/4 for spatially antisymmetric states. Almost all of 
the bound states of HJ are in the Isag ground electronic 
state, which has an attractive potential with a minimum 
energy of E^in = — 2.79eV at an internuclear separation 
R = 1.06 A. The next-lowest electronic state 2pau is repul- 
sive (except for a weak attractive region at large distance 
due to the polarizability of H). 

The radiative association reaction (Eq. |1J begins with 
an H(ls) atom and H^ ion approaching each other. This ini- 
tial electronic state is a superposition of the Isag and 2pau 
states of HJ. Because of dipole selection rules, the system 
can only produce an Hj(lsa"g) ion from the 2pau initial elec- 
tronic state. This state is repulsive, so at low initial energies 
the wave function is confined to large internuclear separa- 
tion R. Wave function overlap considerations then imply 
that radiative association to highly excited Hj(lscrg) states 
is preferred, since these states have significant wave func- 
tions at large R. Conversely, direct radiative association to 
the ground state (A'' = ii = 0) is suppressed. This state of 



HJ can be populated by radiative transitions from excited 
states, but the inversion symmetry of the HJ ion implies 
that these must be electric quadrupole transitions, hence 
they are slow, compared to destruction of Hj. This circum- 
stance results in level populations of HJ that are very far 
from LTE. The non-LTE distribution, with higher-energy 
levels overpopulated relative to the Boltzmann distribution, 
results in photodissociation cross sections that are signifi- 
ca ntly higher tha n the commonly used LTE cross sections 
of lArgvrosI Jl974l) and lstanciL 119941) . 



4.1 Rate equations 

The level populations of H^ are determined by the solution 
of a network of production, destruction, and level-changing 
reactions. Schematically, one may write 

dxi 
'dt 



^ = Si -f ^ RijXj ~ ^ Rj 



(19) 



where st is the rate of production of the ith level of H^ (in 
units of ions per H nucleus per second), Rij is the rate for 
transitions from the ith level to the jth level, and 7i is the 
rate for destruction of HJ ions in the ith level. Note that 
the level index i encodes both the vibrational and rotational 
quantum numbers: i = {v, N). We track all 423 bound levels 
of the ground electronic state, which have quantum numbers 
ranging up to u = 19 and N — 35. 

The source function of Hj comes from the radiative 



association reaction (Eq. 01 and its rate is given by 



(20) 



The level-changing rates include both radiative and coUi- 
sional rates, Rij = Jj'j^'*'' -I- -R'™'^ The radiative rates are 
given by the standard expression 



R 



(rad) 



Ajif{vji)gi/gj 



m < E{j) 
E{j) > E{i) 



(21) 



where Aij is the Einstein coefficient. In principle there is also 
a coUisional term Rlj°^\ which could accelerate H2 produc- 
tion by de-activating H J ions into lower energy levels. These 
lower-energy ions would survive longer since they suffer less 
photodissociation, and hence have a higher probability of 
undergoing charge transfer to produce H2. However, it is not 
possible to produce more H2 molecules by this mechanism 
than there are de-activating H^-H collisions. We will show in 
Section|S|that even if we assume the Langevin rate for charge 
transfer Hj(H,H"'")H2, we find that only ^ 0.005 per cent of 
the H2 is produced via the Hj mechanism. Thus our con- 
clusions about the final H2 abundance are unaffected except 
in the highly unlikely circumstance that the de-activation 
rate coefficient is several orders of magnitude larger than 
the Langevin rate. 

The destruction of Hj proceeds by photodissociation 
(reverse of Eq. charge transfer (Eq. or dissociative 
recombination (Eq.jnj, at a rate 7^, 



Ji ^ Pi + kl'^'^xmn + k^f^'^XeU. 



(22) 



The radiative association and dissociation rates Qi and Pi, 
and the quadrupole Einstein coefficients Aij, are computed 
in Appendix ^ We examine the charge transfer (Section 
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14. 2t and dissociative recombination fSection 14.311 reactions 
in the following sections. 

As described in Section [4.31 we neglect k^'^'^^ since it is 
small during the regime where the mechanism is most 
active; but if it were included, its only possible effect would 
be to further reduce the already negligible H2 yield from this 
mechanism. 

Eq. 1191 possesses a steady-state solution 

x = T-^s, (23) 

where we have written the level populations as a vector and 
the matrix T is given by 



Tij = —Rij + 5ij + ^ Rki^ 



(24) 



which is valid if all the eigenvalues of T are large (fast) 
compared to the Ifubble time. The H2 production rate via 
the mechanism is then 

S+ =na;Hi^fcf'^a;,. (25) 



4.2 Charge transfer 

In order to complete our analysis, we need the rate of the 
charge transfer reaction (Eq. |^ as a function of the matter 
temperature for each level of the initial-state ion. Unfor- 
tunately, ther e are no published computations of the state- 
resolved rates llLepp et al .l2002h . Most pregalactic chemistry 
networks have used the value 



i.(ct) 



(6.4 ± 1.2) X 10 



-10 3 -1 
cm s 



(26) 



measured bv lKarpas et al.l ^1979^ in an ion cyclotron reso- 
nance device; however the depende nce on the tem perature 
and initial state was not determined. iKrsti c (2002) has com- 
puted cross sections for Eq. Q resolved into individual vi- 
brational levels, but was motivated by studies of controlled 
fusion plasmas and so only extends down to thermal energies 
i^ksT^ = 0.1 eV at 2 = 300) for some of the levels. 

Given that the rate for Eq. @ has not been measured 
or calculated accurately for all relevant levels, we have run 
our calculation for three different cases. In case (A), we use 
the lKarpas et all (^7^ rate coefficient ('Ea. l26ll for all lev- 
els; in case (B) we have used the Langevin rate coefficient, 
;.(ct) ^ 2.38 X lO'^'cm^s'^ and in case (C) we have lin- 
early interpolated the values of aVrd from lKrstid (|20o3) and 
obtained a reaction rate by integrating over a Maxwellian 
energy distribution. 



/^(fcsT„)3/2 



T-il/2 -E/kgT^ , CI 

h ' e ' " crUroidi!/, 



(27) 



where E is the kinetic energy of the H-I-Hj system in the 
centre-of-mass frame and tirei is the initial relative velocity of 
H and . For energies less than the lowest tabulated value 
(0.17 eV in the case of the u = state) we have assumed av-ccx 
to be constant. This probably overestimates the reaction 
rate since for most levels (THroi is an increasing function of 
collision energy. 



4.3 Dissociative recombination 

It is possible for dissociative recombination (Eq.|SJ to reduce 
the production of H2. In general dissociative recombination 
can be important if it contributes significantly in Eq. 12211 . 
i.e. if 



,{dr) a:^HI (ct) _ ^-l^{ct) 



(28) 



Of the various models we consider, the lowest value of fc. 



occurs for model A, with fc' 



(ct) 



6.4 X 10"^°cm^s 



(ct) 
i 

At 



z < 600 where the Hj mechanism is most active, we have 
Xe < 10"'' so in order for dissociative recombination to be 
important, k^'^'^^ would have to be 6 x 10~^ cm^ s~^. The 
tabulated rate coefficients for the fir st 78 levels of Ht are 
less than this for 20 ^ r„ ^ 5000 K JSchneider et alJll994l : 
the maximum value in the table is 3.1 x 10~'^cm^s~^ at 
Tm = 20K, and 1.3 x 10"^cm^s"^ for > 100 K); we 
therefore neglect dissociative recombination. This may not 
be a valid approximation at 2 > 600 where Xe > 10"'^, or 
for the highest excited levels of (for which no published 
rates are available). If dissociative recombination is impor- 
tant in these circumstances, the effect would be to decrease 
the II2 abundance at 2 > 300. There would not be a sig- 
nificant effect at lower redshifts, because at 2 < 300 the 
production of H2 is dominated by the H~ and HeH^ mech- 
anisms; is unaffected by dissociative recombination, and 
HeH"'" produces in one of the low-lying states for which 
we have already concluded that dissociative recombination 
is irrelevant. 



4.4 A "two-level" approximation 

Although the above discussion completely specifies the so- 
lution for the channel, it is useful to consider an approx- 
imate "two- level" solution. This serves both as a check of 
the more involved numerical calculations above, as well as 
highlighting the essential physics behind this channel. 

We start by describing the channel by the following 
two level decomposition 

• II2 forms in an excited state Hj(i), at a rate at, 

H + H+^H+(i), (29) 

• The excited ion is then either photodissociated at rate 



A, 

Hj(j) + 7^ H + H^ 



(30) 



or transitions to a LTE distribution via a quadrupole tran- 
sition at rate TZi, 



H+(i)^H^ 



(31) 



Note that although this is a two level description, we still 
consider all energy levels i; they are now however decoupled 
from each other. 

Given Eas. 1291 to 1311 we can write out the differential 
equations describing the time evolution of the abundances 
of these species. Furthermore, as before, Eas. I30l and l31l oc- 
cur significantly faster than the Hubble time, and so, we 
use the steady-state approximation to eliminate any explicit 
mention of the Hj(i) abundance. This yields the following 
equation for the abundance. 
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Table 1. The suppressed association rate as a function of rodshift, 



Ea.l34t. 


z 


logj^Q a' (cm^ s ^) 


z 


logjQ a' (cm^ s ^) 


50 


-21.8824 


300 


-22.1023 


100 


-22.1048 


350 


-22.2281 


150 


-21.9867 


400 


-22.3786 


200 


-21.9676 


450 


-22.5421 


250 


-22.0103 


500 


-22.7089 



da:[Hj 
At 



E 



ailZi 



na;Hix[H+] - /32;[H^ 



(32) 



where /3 is now the LTE photodissociation rate (determined 
by averaging Pi as calculated in Appendix IXlover the Boltz- 
mann distribution of levels at the CMB temperature). Fi- 
nally, there remains the issue of what to assume for TZi ; for 
simplicity, we assume that TZi is the total transition proba- 
bility to lower levels. 



E(j)<E(i) 



(33) 



Note that Ea. l32l now resembles the standard calculation 
with a suppressed radiative association rate. 



EciiTZi 



(34) 



Note that the error in this approximation could go either 
direction. This is because, although our nascent ion is 
likely to be produced in a highly excited state, we do not 
really know without the full multi-level calculation whether 
after it emits its first photon the resulting distribution is 
"more excited" or "less excited" than Boltzmann. We defer 
a comparison of this approximation with the full calculation 
until Section |S| and conclude by tabulating the suppressed 
association rate as a function of redshift (Table 0. 



5 THE HeH+ MECHANISM 

The HeH"*" mechanism is really an additional set of reactions 
that couple to Eqs. Therefore the most straightfor- 

ward way to include it is to add HeH"*" as an additional 
"level" in the network of Section [4. II Due to its large dipole 
moment, HeH"*" has very short-lived e xcited levels (lifetime 
~ IQ-^ s: iRoberge fc Dalgarnolll^82^ . These lifetimes are 
short compared to the photodissociation time from these 
levels or the time between collisions (~ 10* s aX z = 300), so 
the distribution of level populations is determined entirely 
by the radiation field. In the absence of a spectral distor- 
tion, then, it would be permissible to assume that the level 
populations of HeH"*" are in LTE at the CMB temperature 

TcMB ■ 

In reality there is a spectral distortion, however its effect 
on HeH"*" is negligible. The distortion photons that would af- 
fect IIeII+ are those aX E > 1.6[(1 -I- z)/250] eV, where the 
distortion is significant compared to the thermal CMB ra- 
diation; there is roughly 1 distortion photon in this energy 
range per H atom. Since the rate coefficient for Hj" formation 
from HeH+ via Eq. © is 9.1 x 10"^° cm^ s"^ iKarpas et alJ 



[19790, photodissociation of HeH by distortion photons can 
be significant if the photodissociation cross section aver- 
aged over the HeH^ level and distortion photon energy dis- 
tributions is ac ~ 9.1 X 10"^" cm'^ s~^. The cross section 
for photodissociation from the v — Q level of HeH'^ is a l- 
ways ac < 1.5 x 10~^^ cm^ s~ ^ iRoberge fc DalgarndlTgsll . 
The cross section from the higher excited state s v = 7 and 
V = 8 is crc ~ 3 X 10~^ cm^ near threshold iSaha et alJ 
[1978^ and could, in principle, be important, except that 
the fraction of HeH"*" in these states is ^ 1 (they lie 1.56 
and 1.62 eV above the ground state, as compared with 
ksTcMB = 0.07 eV at z = 300). Hence the HeH+ levels 
are in LTE at the CMB temperature Tcmb and HeH^ can 
be followed as an "effective l-level" molecular ion. 

It is straightforward to write down the additional terms 
in Eqs. (I20H22II to take into account the HeH^ contribution. 
The additional source is due to radiative attachment. 



ShcH+ = aHcH+a^pelxHm, 



(35) 



where aHeH+ is the radiative attachment rate coefficient. 
The destruction term is simply the photodissociation rate. 



7hcH+ — /^HcH+ • 



(36) 



Finally the transition matrix R picks up two additional 
terms due to Eq. JSJ. One is the forward reaction term 



Ri 



HoH+ = ^i^2;mn 



(37) 



where ^ is the rate coefficient for Eq. ||SJ and is the 
branching fraction to the ith level of H^. The other is the 
reverse reaction, which can be obtained by the principle of 
detailed balance. 



3/2 



(38) 



R ^ _ R „[B(i)+Oo(HoH+)]/fcBT„ 

2Q(HeH+,r„0 / 

The ratio of reduced masses is MhcH++h/Mho-i-h+ ~ "'Z^' 
and the partition function Q(HeH'*") is obtain ed at the mat- 
ter temperature T™ from lEngel et alj l|2005 [). The bindin g 
energy is _Do(HeH+) = 1.84412 eV iZvgelman et alJll99dl . 
The factor of 2 comes from the ground state degeneracy of 
H (He has degeneracy 1). The remaining additional term in 
the transition matrix is -RhcH+,HcH+ — 0- 

To compute the H2 production, one must also know 

«HoH+ , /3hoH+ , and /H- For aHeH+ > we have 

used the fit bv iGalli fc Fallal Jl998ll to the results of 
iRoberge fc Dalgarnd dl^S^ This rate can in principle 
be increased by stimu lated radiative association. However 
IZvgelman et alJ lll998h found that this increases the HeH"*" 
abundance by < 20 per cent in the redshift range of interest, 
and since our analysis does not change the formation and/or 
destruction mechanisms a similar result would apply to our 
case. We have thus not included a correction for stimulated 
radiative association (we will see that the HeH''' mechanism 
is not a major source of primordial H2, so a correction of 
this magnitude in HeH''" abundance translates into a much 
smaller correction to the net H2 production). For /3hcH+i we 
have used the principle of detailed balance. 



HoH+ — 



Q{ReR+,TcMB) 



/ ^kbTcmb Y 
\ 2Trh^ J 
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xc 



-Do{HcH+)/fcBTcAfi 



oihcH+{Tcmb) 



(39) 



Net H2 production 



here /i is the reduced mass of He and H . The rate c oef- 
ficient fcj^ has been measured by iKarpas et al.l il979l) to 

be 9.1 X 10"^° cm^ s"^ For the f^, which control the dis- 
tribution of rotation-vibration levels in the final state of 
Eq. JHl, there do not appear to be any published measure- 
ments or calculations. We have thus considered two cases: 
one in which Eq. ||HJ populates all energetically available 
levels with their statistical ratios, and one in which the re- 
action always leaves in the ground state of either para- 
H+ (w = 0, = 0; probability 1/4) or ortho-H+ (v = 0, 
A'^ = 1; probability 3/4). While the two cases lead to differ- 
ent level populations of Hj, the total abundances a;[Hj](z) 
and a;[H2](z) are unaffected. This is because, for any of the 
levels accessible to Eq. Q at thermal energies, the most 
probable fate of the ion is to radiate away its vibrational 
energy on a timescale of order 10^ s, and then to undergo 
charge exchange to produce H2 on a timescale of order 10^ s. 
Thus the information about the initial {v, N) distribution of 
the is erased. Photodissociation from the \ow-v levels 
of is strongly suppressed due to lack of wave function 
overlap, and electric quadrupole excitation to high-n levels 
followed by photodissociation is slower than charge transfer; 
thus photodissociation is not effective at depleting the 
produced by Eq. Q. The case where all energetically avail- 
able levels are populated with their statistical ratios will be 
used in the rest of this paper. 



6 RESULTS: H2 ABUNDANCE 



Indiv. mechanisms 
Total Hp (new calc] 
Total (old calc.; 




20 30 50 100 200 300 500 
1+z 

Figure 1. The total abundance of H2 as a function of redshift. 
The thick solid lines shows the new calculation of H2 abundance 
for models A and B (C is indistinguishable from B on the scale 
of the plot). The thick long-dashed line shows the old calcula- 
tion, which overestimated x[H2]. The thin short-dashed lines show 
the individual production mechanisms: H~, H^/HeH"*", and 
(without HeH+). The triangles indicate the results for only 
from the simplified model of Section 14.41 



We now present the results of integrating the production 
rate of H2 from the three major mechanisms. The total pro- 
duction of H2 is obtained by the integral 



[H2]= iS-+S+)dt, 



(40) 



which is shown in Fig. The production rate per Hubble 
time (S- -1-5'+)///' is shown in Fig.|21 The abundances of the 
intermediates H~, Hj, and HeH"*", obtained from Eqs. 1121 
I23II . are shown in Fig. |21 



6.1 H" 

In our calculation, the final H2 abundance is determined es- 
sentially entirely by the H~ reaction sequence. In accordance 
with previous calculations, at high redshift {z > 140) the H~ 
ion is formed mainly through radiative attachment, and de- 
stroyed mainly through photodetachment by thermal CMB 
photons. Note that this is not quite a Saha-type equilibrium 
because of the different matter and radiation temperatures. 
At z ~ 127, the spectral distortion begins to dominate the 
photon spectrum at energies of ~ 1 eV, where the H~ pho- 
todetachment cross section peaks. Below this redshift, H~ 
is still produced mainly by radiative attachment, but the 
destruction mechanism is photodetachment from distortion 
photons. This situation remains until z ~ 67, when most of 
the spectral distortion has redshifted to below the H~ pho- 
todetachment threshold. At this time, the competing pho- 
todetachment and associative detachment (Eq.|5J rates are 



10"' 
10-' 



N 

■.- -fn-6 



CO 



1 10-^° 



10 



11 



Hg production rates per Hubble time 



Indiv. mechanisms 
Total Hp (new calc.) 
Total flj (old calc.) 




(incl. HeH-" 



H2 (without HeH-") 



20 30 



50 



100 
1+z 



200 300 500 



Figure 2. The production rate of H2 via the H~ , H^, and HeH+ 
mechanisms, in units of H2 molecules per hydrogen nucleus per 
Hubble time. The thick solid line shows the new calculation of 
H2 rate for models A, while the thick long-dashed line shows the 
old calculation. The thin short-dashed lines break down the H2 
production into the individual mechanisms. 
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Abundances of intermediate species 




20 30 50 100 200 300 500 

1+z 



Figure 3. The abundances of the intermediates H , , and 
HeH+. The curve is shown for both models A and B. Note 
that a;[H~] peaks at ^ ^ 62 when H2 production is also maxi- 
mized. 



similar (~ 10"^° s~^). This is also the era of peak produc- 
tion of H2: at higher redshifts the H~ ions are destroyed 
before they can react with H to produce H2, while at lower 
redshifts less H~ is produced due to the lower density of 
the universe and the decrease in reaction rate at low r,„. 
Ultimately we find that the total amount of H2 produced 
via the H~ mechanism, J S-dt, is only 6 x 10~^ instead of 
2.2 X 10~^ as found in the standard calculation. 

As noted in Section |3 there is some uncertainty in the 
rate for associative detachment, Eq. Glover ct aj ll2006t) 
argued that the rate coefficient ^2 could plausibly be varied 
between 6.5 x 10"^° and 5 x 10"^ cm s . We have re-run 
our analysis using these values and find that the final H2 
abundance varies from x[K2] = 4.1 x 10~^ for the lowest 
value of k2 to 9.4 x 10~^ for the highest value. Thus our 
fiducial estimate of 6 x 10~^ should be considered uncertain 
by a factor of ~ 1.5 in either direction. 

There is also a large uncertainty in the mutual neu- 
tralizat ion rate, with some results (e.g. the experimental 
work of lMoselev et aLHigTOf) being up to an order of magni- 
tude higher than the fits used here. Using the higher mutual 
neutra lization rate determined by the fit of lMoselev et all 
il970l) . we find that the final H2 abundance decreases from 
6.0 X 10~^ to 5.7 X 10~^. Thus the uncertainty in the mutual 
neutralization rate does not have a significant effect on the 
final H2 abundance. 



6.2 H+ and HeH+ 

We find that the reactions involving positive ions (Eqs.UHHJ 
do not contribute significantly to the final H2 abundance. 
Nevertheless these reactions operate at earlier times than 
H~ due to the higher binding energy of HeH"*" and com- 
pared to H~, and so they dominate the H2 production at 
z > 144. The total amount of H2 produced by these reactions 



is J S+dt = 5 X 10 ®, with almost all of this contributed by 
HeH+. 

Our computed HeH'*' abundance (see Fi g. are very 
similar to those obtained in previous works ilGalli fc Fallal 
Il998l: IStancil et aLllToil) . In contrast, we find much lower 

ion abundances (and H2 production rates via Hj) than 
in the standard calculation. This is a consequence of the non- 
Boltzmann level populations in HJ. An example of these 
level populations for model A at z = 300 is shown in Fig. |31 
(Results are qualitatively similar for model B, with the main 
quantitative difference being that the lowest-energy states 
are less populated.) As shown in the figure, the highest- 
lying levels are nearly in Boltzmann equilibrium since the 
reaction 

H + H+ ^ Hj(lscr9,i;iV) -f 7 (41) 

is fast compared with the electric quadrupole transitions 
in HJ. (The slight deviation from Boltzmann equilibrium 
among the high-lying levels in the figure is due to the dif- 
ference between matter and radiation temperatures.) The 
populations of the lower-lying levels are determined by a 
combination of quadrupole radiative cascade rates, sourcing 
by HeH^ through Eq. ©, and (at the lowest energies) some 
quadrupole excitation by the CMB. Note that even in the 
ground vibrational state {v — 0), the rotational levels of 
never come to thermal equilibrium with the CMB because 
the radiative rates (which may be as long as ~ 10^" s) are 
slower than the timescale for charge transfer (about 3 x 10* s 
for model A at z = 300). 

The overall production of H2 via the HeH^/Hj mech- 
anisms IS found to be x[H2] = 4.8 x 10"® at z ^ 20. Of this, 
only a small fraction (a:;[H2] = 2.6 x 10~^^) is formed if we 
artificially turn off HeH"*". (These numbers are for model A, 
models B and C give numbers that are < 10 per cent higher.) 
The direct production of H J via Eq. (01 only contributes sig- 
nificantly at z > 300 where it is possible to directly produce 
tightly bound states by radiative association. 

The two-level model of Section 14.41 does quite well at 
reproducing the final H2 abundance due to (it yields 
2.2 X 10"" versus 2.6 x 10"" for the full multi-level calcu- 
lation). It is not quite so good at reproducing the redshift 
history: at very high redshifts where can be excited out 
of the ground state and photodissociated, it underestimates 
H2 production because the photodissociation rate of the 
ion after it has reached a low-lying level is actually less than 
the LTE rate. 



7 HEATING 

The principal reason for interest in the molecular abun- 
dances in primordial gas is the possible effect on the heating 
and cooling rates. Before the formation of collapsed struc- 
tures, the gas is colder than the CMB and hence the effect 
of molecular transitions is to heat the gas; after collapse the 
gas is compressed and shock-heated to above the CMB tem- 
perature, at which point molecules cool the gas. The latter 
regime is dominated by H2 molecules formed in the collapsed 
gas clouds themselves rather than primordial molecules, and 
hence it is not considered here. The molecular heating in 
the pre-coUapse regi me was first inv e stigate d by|Pm_£t_aL 
119931 . and later bv iHaiman et all (|l^63), ICaUi fc Fa~ 
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H2 level populations at z=300 
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Figure 4. The level populations in at 2 = 300 for model 
A. The horizontal axis shows the level energy with E = de- 
fined as the dissociation limit into H+H+. The actual populations 
for each of the 423 rotational- vibrational levels are given by the 
crosses. The Boltzmann equilibrium populations oc e~^^^^^ at 
the radiation and matter temperatures are shown by the straight 
lines. Note that the levels close to the continuum are approxi- 
mately in equilibrium, but the low-lying levels are highly under- 
populated relative to Boltzmann. The levels more tightly bound 
than E{vN) < — 1.844 eV have higher populations because they 
can be formed from HeH"*" via Eq. JSJ. At this redshift the ra- 
diation and matter temperatures are kgTcMB = 0.0708 eV and 
^BTm. = 0.0662 eV, respectively. 



HHi), and iFlower fc Pineau des ForetsI (|200(t . The pur- 
pose of this section is to revisit the heating rate, taking 
into account the revised molecular production rates, and in- 
cluding the release of heat in c hemical rea ctions. The chem- 
ical heating was considered bv lPuv et alJ (il993.) and found 
to be "negligible," but no quantitative result was given; in 
any case we have found large changes in some of the rele- 
vant reaction rates and it is only prudent to re-evaluate the 
chemical heating term. 

The overall temperature balance equation is 



dTrn 

dt 



-In Ira H (J-CMB ~ J- rr 



3kBXt 



{Qff + Orec + 



^^,(42) 
xt at 



where the first term represents the adiabatic expansion, the 
second term represents Compton heating, the third term 
represents heating by molecular rotation and chemical re- 
actions, and the fourth term takes account of the changing 
number of translational deg rees of freedom among which 
kinetic energy is distributed llPuv et alJll993ft . In this equa- 
tion, ot is the Thomson cross section; abb is the blackbody 
radiation coefficient; me is the electron mass; xt is the num- 
ber of gas particles per H nucleus; Q// is the free-free heat- 
ing rate; Qrec is the net heating due to recombinations (note 
that this is negative); Qrot is the heating by molecular rota- 



tion per H nucleus per unit time; and Qc is the heating by 
chemical reactions per H nucleus per unit time. 



7.1 Rotational transitions 

In this section we evaluate the heating in the rotational 
lines of H2. It is found to be negligible, and is much smaller 
than the direct heating produced by the chemical reactions 
that generate II2. Determining the rotational line heating 
requires determining first the level populations and then 
computing the heating rates by considering all coUisional 
excitations and de-excitations via H and H+. 

The H2 molecule has two sets of levels with different 
nuclear spin properties, namely para-H2 (even angular mo- 
mentum J, and total nuclear spin 7 = 0) and ortho-Il2 (odd 
J; I — 1)} In the high-redshift universe, the timescale for 
the electric quadrupole radiative transitions that connect 
H2 levels with the same nuclear spin is fast compared to 
the time scale for collis i ons. F or example, using the radiative 
rates of iTurner et al.l il977l) and the coUisional rates de- 
scribed below, we find that the collisional-to-radiative tran- 
sition rate ratio at z = 250 is 0.04 for J — 0; this rate 
is even less for higher rotational levels or lower redshifts. 
Therefore we may treat the H2 molecule as an effective two- 
level system, tracking separately the abundances of para-H2 
and ortho-H2. Within each set of levels (even or odd J), the 
populations are assumed to rapidly thermalize to the CMB 
temperature, a;[Il2, J] oc gje~^''^^^'^'^'''^ . Here the degen- 
eracy is gj — (2 J + 1) /4 for even J or gj = 3(2J + 1) /4 for 
odd J. The evolution equation for the fraction of the H2 
molecules in the ortho form is 

dt 



(•^ prod 



.^)_lna;[H2]+7p^o(l- 



- (43) 



where .7-"prod is the fraction of H2 that is produced in the or- 
tho levels; ^p-,o is the para-to-ortho transition rate (in s~^); 
and 7o^p is the ortho-to-para transition rate. The transition 
rates are 

7p_o = n (x[H+]{ap^oi;)jj2+H+ + xm{(Tp~.ov)ii2+ii) , (44) 

where the averages are taken over both thermal velocity at 
temperature Tm and the Boltzmann distribution of H2 ro- 
tational levels. A similar equation holds for 7o^p. We have 
taken the cross sections for H2-| -H^ from lG erlichI lll990ll . For 
H24-H, we have used the fits by 'Flower* ('19971), with the cor- 
rection for reactive scattering by Lc Bou rlot et al.l il999l) . 
Also, since H2 is usually produced in a highly excited rota - 
tional state when it forms from H~-|-H iLaunav et al .119911) . 
we have assumed here that J^p^od = 3/4 as would be sug- 
gested by the nuclear spin degeneracy. Once the level pop- 
ulations are established, the heating rate is computed by 
assuming that an amount of energy Ej — Eji is added to 
the translational degrees of freedom of the gas for each col- 
lisional J J' transition: 



Qrot 



E 

.J.J' 



x[R2,J\{Ej - Ej,)n[x[R-^\{aj_j,v)^^+^+ 
+3::hi(o-./^./''!;)h2+h). (45) 



^ The orbital angular momentum N and the orbital plus elec- 
tron spin angular momentum J are identical for the H2 electronic 
ground state since the term symbol is '^S^. 
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Our treatment of H2 cuts of f the rotationa l levels at Jmax = 
9, the highest level for which lGerlichl Jl990l) provides H2+H+ 
cross sections; but we have found that at z < 400 there is 
a < 10 per cent change in the heating rate if we cut off the 
levels at Jmax = 7 instead. 



7.2 Chemical heating 

We next consider the chemical heating from each of Eqs. Q- 
The first contribution comes from the formation and de- 
struction of H[^, Eq. Q. The rate of loss of kinetic energy 
from the forward reaction is kiE{Tm)xeXuin, where E{Tm) 
is the mean energy of the incident electron participating in 
the reaction. This mean energy is given by multiplying the 
Maxwell distribution against the radiative association cross 
section, 



Qc2 = (0.91eV)A;2a;[H"]a;Hin. 



(48) 



E = 



-E/kBTm 



AE 



Eara{E)e-'^/^BT,^AE 
A 



knT, 



'ATm 

A 

'ATm 

' (2 



■In 



EGra{E)e 

ln[T^/'fci(T™)] 
din fci 



-E/kBT„ 



AE 



AlnT„, 



(46) 



where in the last line we have re-expressed the energy inte- 
gral in terms of the thermally averaged rate coefficient. By 
the principle of detailed balance, the average energy of the 
photodetached electrons from the thermal part of the reverse 
reaction is given by the same function, E{Tcmb)- For the 
electrons that are photodetached by the spectral distortion, 
the kinetic energy input is hv — i3(H^). Thus the overall 
heating term associated with Eq. Q is 



Qci 



3 din fci(r^ 
2 ^ A\nTm 



kBTmkix^xuin 



+ 



+ 



^[H-]{ 



3 A\nki{TcMB] 
2 d In TcMB 



kBTcMBk_-^ 



(th) 



B(H- )/h 



r{u)a^i{u)[hiy-B{R-)]^y (47) 



The second contribution is the heating due to as socia- 
tive detachment, Eq. Q. As noted bv lPxiv et alJ 1^23)) this 
reaction is exothermic by 3.72 eV; however it is not correct 
to set the heating term from Eq. (|5J equal to 3.72 eV times 
the reaction rate, because this energy yield is distributed 
among both the kinetic energy of the ejected electron (which 
couples to the gas temperature) and the rotational and vi- 
brational degrees of freedom of the H2 molecule (which are 
radiated away on timescales of ~ 10® s, i.e. much shorter 
than the collision timescale). In order to compute the heat- 
ing term, one needs to know the mean excitation energy of 
the final-state H2 molecule. This can be determined from the 
(f, J)-resolved cross sections for associ ative detachment of 
H and H~, which were computed bv iBieniek fc Dalgarnd 
il979l) at an initial state energy of 0.0129 eV (which is 
roughly the initial energy of interest); this yields a mean 
excitation energy of 2.81 eV. By subtraction we estimate 
that 0.91 eV of energy is available to heat the gas. Thus we 
set 



A third contribution to the heating rate is from mu- 
tual neutralization, Eq. Q. The branching fraction for this 
rate is small, but it is strongly exothermic (12.84 eV) and 
thus provides a potentially large amount of energy. The fi- 
nal state products are If atoms, which have no rotational or 
vibrational degrees of freedom. They do however have elec- 
tronic degrees of freedom: it is energetically possible for the 
final state to be H(ls)-|-H(n/) for any ^ I < n ^ 4. At 
low energies, nearly all of the neutralizations go to n = 3 
( Fussen fc Kubach,1986.) . so the energy released into trans- 
lational degrees of freedom is 0.76 eV. The excited H atom 
decays by emitting Ha, Lyman, and/or 2-photon continuum 
radiation, none of which can heat the gas. Therefore the mu- 
tual neutralization heating term is 



Qc3 = (0.76eV)fc32:[H-]a;[H+]n. 



(49) 



We also consider the heating and cooling from the Hj 
and HeH"'' reactions. The contribution from Eq. I^J is ob- 
tained by inserting the translational energy E — E{vN)+hv 
in the integrals for radiative association (Eq. lAlOt and 
photodissociation (Eq. IA5II . The heating from the charge 
transfer reaction, Eq. (|^, depends on the final rotational- 
vibrational state of the H2 molecule, which is not resolved 
in our code. We have assumed that the H2 is produced in 
the ground state, which maximizes the heating from this 
reaction since it implies that the entire energy yield of the 
reaction is available to heat the gas; however Eq. JSJ is not a 
significant source of chemical energy anyway. The radiative 
association and dissociation of HeH"*" (Eq. |7| was treated 
using the obvious analogue of Eq. 14711 . without the spec- 
tral distortion since as argued earlier it is unimportant for 
HeH"''. The contribution to the heating rate from all of the 
HJ and HeH^ reactions is negligible. 



7.3 Results 

The results are shown in Fig. |5] it is seen that the chemical 
heating is very small, of the order of 10~*. The principal 
source of heating is the H~ sequence of reactions, which 
provide a peak in the heating at 2: ~ 120 from the pho- 
todetachment process, and a second peak at z ~ 70 due to 
associative detachment. There is even less heating at earlier 
times, and in any case heating before Compton freeze-out 
[z ~ 200) will be erased. 

We find that the heating due to rotational transitions in 
H2 is negligible, in qualitative agreement with the results of 
[Plowcr fc Pincau dcs Forcts (2000). The fractional heating 
rate peaks at ~ 3 x 10~® at z « 50, shortly after the peak 
in H2 production rate (there is less H2 at earlier times, and 
the decreasing density and temperature suppress coUisional 
rates at la ter times). The d ifference from some previous re- 
sults (e.g. iPuv et"al]ll993D is due partially to our reduced 
H2 abundance but also due to reduced J-changing collision 
cross sections. In particular, collisions with H""" dominate 
the transfer of rotational energy into translational degrees 
of freedom. While there remains considerable uncertainty in 
the potent ial for the H2-I-H system as noted by several au- 
thors (e.g. iGalh fc Palialll998l : iFlower fc Pineau des ForetsI 
l2000h . we find that the cross sections would have to be in- 
creased by a factor of 40 to make rotational heating con- 
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Figure 5. The heating rate due to the H~ and H^/HeH+ se- 
quences of reactions, and the H2 rotational transitions, for Model 
A. What is shown is the heating input per Hubble time, Q/H, 
divided by the thermal energy of the gas, 3kBTmXt/2. The H~ 
reactions affect the gas temperature at the ~ lO"**, and re- 
actions have even less influence. The H2 rotational line heating is 
always small compared to the chemical reactions. 

tribute at the 10"^ level (2Q^/3//fcsT„ = 10"^), and by a 
factor of 500 to make it contribute at the 10~^ level. Thus 
the conclusion that rotational heating is negligible is rel- 
atively insensitive to the remaining uncertainty in H2+H 
cross sections. 

Chemical and H2 rotational heating are clearly very 
small effects and the only proposed observational technique 
that could achieve the 10"* level of accuracy, even in p rin- 
ciple, is 21-cm tomography jLoeb fc Zaldarriagall20(M) . In 
practice the experimental challenges of detec ting the cosmo- 
logica l 21-cm signal are substantial (see e.g. iBowman et alJ 
I2OO6I for a recent review) , and observing the pre-reionization 
epoch to such high accuracy should be considered a goal for 
the distant future. Until such data are available the chemical 
and rotational heating of the pregalactic gas can be safely 
neglected. 



(via the H~ sequence, Eqs.^l^J rather than rotational lines; 
however the effect is only of the order of 10~*. In princi- 
ple the proposed 21-cm tomography of the pre-reionization 
Universe could reach the sensitivity at which primordial H2 
becomes important, since it is sensitive to the temperature 
of the gas and has many more than (10~'')~^ ~ 10* modes. 
However when assessing the prospects, it should be remem- 
bered that the high-redshift 21-cm signal has not yet been 
detected, and measurements at the 10~* level are clearly 
very far in the future. 

Aside from H2, there are other molecules with rota- 
tional lines such as HD and LiH, which could conceivably 
have been formed in the early Universe and played a role in 
the thermal balance. The treatment of these trace molecules 
is beyond the scope of this paper, but HD in particular 
may warrant further study as it has been found t o be a sig- 
nificant heating source in some past works (e.g. IPuv et al] 
[1993). Since the main route of formation of HD is via the 
reaction H2(D+,H+)HD jGaUi fc Pallalll99i) . any analysis 
of HD must incorporate the revised H2 calculation presented 
here. 

Finally, one could ask whether the H2 suppression mech- 
anisms discussed here - the spectral distortion and non- 
equilibrium populations in Hj - have a significant effect 
on H2 cooling of protogalaxies. In the case of the spectral 
distortion, we have seen that at mean density the photode- 
tachment of H~ becomes unimportant at low redshift, since 
H~ ions undergo a chemical reaction (usually associative 
detachment) before being destroyed by radiation. For ex- 
ample, the branching fraction for H~ at mean density to be 
destroyed by radiation is 0.13 ai z — 40, 0.09 at z = 30, 
and 0.05 at 2 = 20. In overdense gas clouds at 2 < 40 the 
coUisional reactions are faster and we conclude that the spec- 
tral distortion should be negligible. We have tried running 
our HJ code for overdense conditions at low redshift (e.g. 
z = 20, 5i, = 10*, T = 10^ K) and find that the H+ levels are 
still far out of equilibrium, with the lowest levels underpop- 
ulated relative to the Boltzmann distribution (at Tcmb) by 
many orders of magnitude. However in such clouds H~ is 
likely to bo a more important source of molecules than Hj 
(see e.g. [Teg mark ot al. 1997). Therefore we do not expect 
our changes in H^ physics to have large consequences for 
the cooling of the first collapsed objects in the universe. 



8 DISCUSSION 

We have reconsidered the production of H2 molecules in the 
pregalactic medium. In contrast to previous studies, we have 
included the spectral distortion in our analysis, and resolved 
all 423 rotational-vibrational levels of the HJ ion. We find 
that in the level-resolved analysis, the Hj reaction pathway 
is greatly suppressed because newly formed H J ions are pho- 
todissociated before they can decay to the ground state or 
undergo charge transfer to become H2 molecules. We also 
find that the H~ ion is easily destroyed by spectral distor- 
tion photons at 2 > 70, so that the production of H2 by this 
pathway is suppressed relative to the standard calculation. 
We obtain a final H2 abundance a;[H2] = 6 x 10~^ assuming 
standard cosmology. 

Unfortunately, the primordial H2 molecules will be very 
difficult to detect. The main effect of H2 on the thermal his- 
tory of the gas actually comes from the formation process 
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APPENDIX A: THE HJ MOLECULAR ION 

The bound states of HJ(1s(T5) are solved within the Born- 
Oppenheimer approximation, which for a diatomic molecu- 
lar ion gives 

*(r,fl)=x(r|i?)5^yiVM„(A), (Al) 

where r is the electron position, R is the internuclear sep- 
aration vector, is the Iso-g electronic wavefunction, 
<1>(_R) is the radial wavefunction, and R is the unit vector in 
the direction of R. The radial wavefunction is a solution to 
the Schrodinger eigenvalue equation. 



V{R) + 



h'^N{N + l) 
2R? 



^{R) = E${R),{A2) 



where E is the bound state energy and ^ is the re- 
duced mass of HJ. The potent ial V(R) is i nterp olated 
from the tabu l ated values of iBates et alJ lll953f) and 
iMadsen & Peekl il97lfl . The radiative transitions among the 
levels of Hitflsgg) ar e treated according to the method of 
iBates & Poot j lll953l) . who give the spontaneous rates Aji 
in terms of the wave functions^ and the quadrupole moment 
M{R) of the H J ion. The q uadrupole momen t is interpolated 
from the computations of iPeek et alj il979l) for R 12ao, 
and using the asymptotic formula M{R) ~ R^ / A — ada%R~^ 
for R > 12ao (where ad ~ 9/2 is the dipole polarizability of 
H). 

The rate for the radiative association reaction (Eq. |^ 
is required for our calculation. The individual rates to and 
from each rotational-vibrational level of HJ as a function 
of matter and radiation temperature are not tabulated in 
the l iterature, but can be derived using standard me thods 
(e.g. IZveelman & Dakarnol [l990l: IStancil et ai]ll993D . For 
the dipole radiative dissociation, one solves for the unbound 
final state Eq. <A21 using the 2pff„ potential and the energ y 
normalization (see e.g. Eq. 4 of lZveelman & DalgarnJl990ri . 
Then the radiative dissociation cross s ection from the vN 
bound level with energy E{vN) < is (lDunrJll968h 



2 Note that P{R) in iBates fc Pootj ^953) is equal to our 
■S>{R)/R. 
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3c 



(A3) 



where hu is the photon energy and 



M2 = (2iV + l)-^ 
+(iV + l) 



iV / <i>{R)D{R)fN-i{R)dR 



^iR)D{R)fN+i{R)<iR 



(A4) 



Here D{R) is the dipole matrix element between Isag and 
2pau states, and /jv±i(-R) are the unbound wave func- 
tions that solve Eq. IIA2II with rotational quantum num- 
bers A'^ ± 1 and energy E{vN) + hv. (Only unbound states 
where N' = N ± 1 contribute due to dipol e selection ru les.) 
The dipole mome n ts are interpolated from lBate^ il95 and 
iRamaker fc Peek! Jl97 3^ for R < 20ao, and the large-i? 
asymptotic expression (Ramakcr & Peck 1972) is used for 
R > 20ao. 

The net radiative dissociation rate, including both ther- 
mal and nonthermal contributions to the radiation field, is 
then given by integrating over all photon frequencies, 



n - 8^ 

PvN — 



E{vN)/h 



ahu/kBTcMB — 



(A5) 



-E(vN)/h 

The spontaneous radiative association cross section is 
related to the radiative dissociation cross section by phase 
space factors. 



roH-t-H+ 



ard,vN[E - E{vN)], 



(A6) 



where w represents the number of quantum states per unit 
wavenumber k per unit volume for either Hj(t)A) + 7 or 
H + . These are given by 



^H+(i,iV)+7 

and 



: 4{2N + l)g'^. 



A-rvk 



■uj 



H+H+ 



H+H+ 



(27r)3 



(A7) 



(A8) 



where fc^, = 2-kv jc is the photon wavenumber and = 
yJ'lyi.E jh is the relative wavenumber of H and H"^. In ad- 
dition to the usual factor 47rA;'^/(27r)'^, Eq. IIA7II contains a 
factor of 2 from the electron spin, a factor of 2 from the 
photon polarization, and the nuclear spin-degeneracy factor 
Snuci Eq. (IA8II contains a factor of 2 from the electron spin. 
The net radiative association rate coefficient to the vN level 
is then given by the usual integral of cross section times 
velocity over the Maxwell-Boltzmann distribution, 

2^/2 



Ear 



iiE)e- 



-E/ksTrr 



[l + f{u)]AE. (A9) 



Here f{u) is the photon phase space density at frequency 
u = [E — E{vN)]/h, and it has been included to take into 
account stimulated radiative association. The phase space 
density from spectral distortion photons is ^ 1 and can be 



neglected, so we appropximate f{v) by a blackbody. Sub- 
stituting Eqs. IIA6HA8II into Eq. <A9t then gives the final 
result, 



OivN 



2V2(2Ar+l)ff;,,fe3 

7rl/2c2(^ifcflT™)3/2 



^ardiv) 



^[hv+E(vN)\/kBT^ 
1 — Q—hf/kgTQ^.js 



Au. (AlO) 



-E{vN)/h 

In this equation g'^^^ — 1/4 for even A*" and 3/4 for odd N. 
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